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Abstract 

We briefly review the oscillator and spin bath models of quantum environments, 
which can be used to describe the low-energy dynamics of open quantum systems. 
We then use them to discuss both the mechanisms causing decoherence at low T, 
and the dynamics of this decoherence. This is done first for a central 2-level sys- 
tem coupled to these environments- the results can be applied to the dynamics 
of quantum nanomagnets and superconducting SQUIDs, where large-scale tunnel- 
ing of magnetisation and flux take place. Decoherence in these systems is caused 
principally by coupling to electrons and nuclear spins- the spin bath couplings are 
particularly dangerous at low T. 

We may also generalise these models to discuss quantum measurements in which 
the measured system, the measuring apparatus, and the environment are treated 
quantum mechanically. The results can be used to calculate the dynamics of coupled 
SQUIDs and/or nanomagnets, in which one acts as a measuring apparatus and the 
other exhibits large-scale quantum superpositions. The same model can be used to 
describe coupled qubits. 
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1 Introduction 



The physics of "open" quantum systems encompasses a vast array of phenom- 
ena, ranging from sub-nuclear to cosmological scales. Nevertheless, in this the- 
matic kaleidoscope there are a few central problems, whose generality makes 
them relevant at almost any scale. Here we look at one such problem - fo- 
cussing on low energy scales (ie., low-T) where clear experimental tests of the- 
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ory arc possible. The problem is that of decoherence - where it comes from, 
and how to understand its effect on the dynamics of quantum effects involving 
many degrees of freedom, at both microscopic and macroscopic scales. 

Discussions of decoherence have often been associated with rather ineffable 
questions about quantum measurements. quantTim cosmology, and quantum 
computation and some of this literature has tended to wander around in circles, 
saying little that might be quantitatively testable in the real world. This has 
given the subject an eccentric reputation, which is unfortunate, since many ex- 
periments are now searching for (and in some cases finding) macroscopic quan- 
tum phenomena in magnets and superconductors. At the same time entangle- 
ment experiments in optical cavities are showing multi-atom entangled states. 
In both kinds of experiment the crucial stumbling block preventing progress 
is decoherence coming from, amongst other things, coupling to environmental 
modes. The decoherence problem is thus an "open quantum systems problem" 
par excellence; it provides an extremely severe test of our understanding of the 
quantum dynamics of complex systems. Since the same physical mechanisms 
responsible for decoherence are involved in many other processes of relaxation 
in the quantum regime (as well as in the crossover to the classical regime), a 
lot turns on these developments. 

In what follows, we first explain the decoherence problem, and describe some 
models used both for this problem and for more general discussions of quantum 
relaxation. We also describe the connection to experiments, particularly those 
in magnetic and superconducting systems. In both cases we explain why the 
most severe source of decoherence in the low-T limit will often be the "spin 
bath" environment of localised modes (defects, spin impurities, nuclear spins) 
coupling to the collective coordinates (magnetisation, flux) of interest. 

It is then but a small step to a discussion of the quantum measurement prob- 
lem, for the simple reason that one can also use superconducting SQUIDs and 
nanomagnets as measuring systems themselves. Thus we also grasp the nettle 
and consider models in which the entire set-up (quantum system, apparatus, 
and environment) is quantised. This not only allows us to calculate the quan- 
tum dynamics of such a set-up, but also to discuss, in a more quantitative way 
than usual, the measurement problem itself. These same models also describe 
a pair of coupled mesoscopic qubits- and we discuss why decoherence typically 
has a much worse effect on the mutual coherence between these than it would 
on each qubit by itself. 

Decoherence Puzzles: The basic idea of decoherence as a solution to the 

measurement problem has been around for some 50 years [1,2]; and as a phys- 
ical phenomenon it has been noted since before Quantum Mechanics. At first 
glance it seems almost trivial- most physicists and chemists are familiar with 
the way in which the phase interference in a "2-slit" kind of experiment is 
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destroyed. So, why care about decoherence now? The short answer is that we 
now need to know rather urgently what mechanisms cause decoherence, par- 
ticularly at the mesoscopic and macroscopic scales, or in entangled systems 
like quantum computers- and we find that there are many basic things that 
we don't understand. One way to introduce the problems is in the form of a 
set of "decoherence puzzles" or questions, which are clearly down-to earth and 
within the purview of physicists and chemists. Here are a few (some of intense 
current interest): 

(1) Mechanisms of Decoherence : What are the mechanisms of decoherence 
in nature and how do they relate to those governing quantum relaxation ? 
This problem is very difficult because decoherence is extraordinarily sensitive 
to the different environmental couplings in nature. Thus ad hoc theoretical 
models, devised for their solvability rather than their realism of generality, 
will not do here. We are entering a particularly exciting time for this puzzle, 
since experiments arc now directly addressing coherence phenomena at the 
mesoscopic and macroscopic levels. 

(2) Measurement Puzzle : Does decoherence "solve" the measurement problem 
(as many [1-3] have advocated)? How then must we understand measurement 
operations in systems where macroscopic coherence exists ? 

(3) Dynamics of Decoherence : How does decoherence (in the relevant compo- 
nents of the relevant density matrix) evolve in time? How does decoherence 
affect N-particle entanglement (eg, in an A^-qubit quantum computeri [4]) and 
how do N-particle decoherence rates depend on N 7 How can the dynamics 
be controlled/tested experimentally ? 

(4) The T(i) puzzle: Why, in defiance of very basic theory [5] docs the decoher- 
ence rate r^^(T) for electron propagation in mesoscopic conductors appear to 
saturate [6] to a finite constant at low T ? Does a similar saturation exists in 
superconductors or magnets ? How can t^^{T — > 0) be controlled ? 

(5) Quantum- Classical Crossover. This is induced by raising cither tempera- 
ture, or external fields, or temperature. The environment inducing the tran- 
sition is usually neither in internal equilibrium nor in equilibrium with the 
system making the crossover. Yet this departure from equilibrium is actually 
crucial (in, for example, many discussions of quantum measurements); how 
can it be treated properly? 

(6) High Energies: What is the role of decoherence in the early universe [7] 
(another possible quantum-classical crossover), or in black hole physics [8] ? 
Can high-energy sources of decoherence up to or beyond the Plank scale have 
an effect at low energies [9]? How can this be tested experimentally ? 

Note that these puzzles are all more or less related, and that they span many 
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different areas in pliysics (witli obvious repercussions in pfiysical cfiemistry) . It 
is also clear that the first puzzle, about the mechanisms governing decoherence, 
is quite central, and is the sort of bread-and butter question that physicists 
and chemists deal with every day. Perhaps the first thing to understand is 
why it is possible to approach this puzzle in a unified way (rather than by 
piecemeal discussion of many different systems), and so we start with this. 



2 Spin Baths and Oscillator Baths 

Theorists working on the quantum dynamics of open systems like to have 
models that are both realistic (realistic enough to survive experimental tests 
on particular physical systems) and also general (so that the results are generic, 
or at least apply to a large class of physical systems). One way of deriving 
such such models is to divide the quantum environment into extended and 
localised modes, and then proceed as follows: 

(i) Extended Modes: In a box of volume V containing Ng ~ 0{V) extended 
modes below an ultraviolet cutoff Ec, the coupling between these modes and 
some "central system" (collective coordinate, etc) of interest is necessarily 
~ 0{N~^^'^). Then for large No one can treat perturbatively and map to an 
environment of oscillators [10,11]. The general model is then described by 
a Hamiltonian H = Ho{P,Q) + H°^^ + Hos{{:>^q,Pq}) where Ho{P,Q) is the 
Hamiltonian of the "central system" we are interested in, having canonical 
coordinates P, Q, and the oscillator terms are [10,12] 

1 No 

Hos^^Y^^qi-kl+U^y,) (1) 
^ 9=1 

No 

HZt-T.[FM^'i + G,{P)p,] (2) 

q=l 

where the {x^, p^} arc canonical coordinates for the extended environmental 
modes (with q = 1,2, ... No). The linear couplings Fg,Gq ~ 0{N^^^'^), so 
that the bath modes have individual effects ~ 0{l/No) on the central system 
- summation over these gives a well defined thermodynamic limit independent 
of No as No oo. 

The derivation of oscillator bath models generally uses some kind of Born/Op- 
penheimer-like assumption that one can separate "slow" system modes (P, Q) 
and "fast" environmental modes {xy,Py}; such derivations are of course very 
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familiar to quantum chemists. To be valid they require that 



K,,'(^,Q)<K-a;,0 (3) 

(or an appropriate generalisation to finite temperature T) for all q, q', P, Q 
of interest, where the matrix element at a given P, Q is given, in terms of the 
full multi-dimensional coordinate space X of the environment, by 



(p,Q) = / dx<f>;{x)HZM^) (4) 



where the wave-functions ^^(X) are the eigenstates of the environmental 
Hamiltonian (for more details sec, eg., Caldeira and Leggett [11]). Note that 
this restriction is not nearly as bad as it seems since one has considerable 
freedom in the choice of the environmental modes- indeed the {cUg}, as well 
as the couplings in (2), typically depend on T (because, eg., one chooses a 
T-dependent effective Hamiltonian for the environment). 

However, there are obvious cases where the assumption (3) breaks down. In the 
case of extended environmental modes a well-known example is the Quantum 
Hall system, where the environmental Landau levels are all degenerate- and the 
effect of H°^^ is typically to break this degeneracy [13]. However the breakdown 
is much more common in the case of localised environmental modes, to which 
we now turn. 

(ii) Localised Modes: In the same volume V there will also be be many 

localised environmental modes, associated with defects, impurities, param- 
agnetic spins and the nuclear spin bath. Each of these has a finite Hilbert 
space in the energy range < Ec, and can thus be mapped to a spin degree of 
freedom, very often a spin— 1/2, or two-level system. Some set {(Tk} of these 
with k — 1,2, . . . Ng, will also couple to the central system, giving extra terms 
H-^^ + Hgp in the total Hamiltonian, where [14] 



= 5: h, • + ^ V^,^,atal, (5) 

h h h' 

HZ = E [fHQH' + (FkiP, QW' + h.c.)] (6) 



The important point here is that there is no requirement that the couplings 
{fJ,F^^} be 0{N-^/^y, indeed very often they are independent of A^s ! A 
simple example is the hypcrfinc coupling, crucial in magnetic systems since it 
couples the macroscopic magnetisation to the nuclear spins- this is obviously 
independent of Ng. Note that the couplings are not necessarily small either; 
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for example, hyperfine couplings can be as large as 0.5K each, and individual 
exchange couplings to paramagnetic impurities can be much larger again. Now 
2 crucial points arise, viz., 

(i) The intrinsic dynamics of the environmental modes, coming from Hsp in 
(5) is now governed by energy scales |hjfc|, \V^^\ which are often much weaker 
than the couplings Fi^}, putting us then in an "anti- Born Oppenheimcr" 
limit, where the environmental modes are strongly slaved to the system co- 
ordinates P, Q. This is because the modes, having localised wave-functions, 
usually couple very weakly to each other (via, eg., residual dipolar couplings). 
We may again take magnetic systems as an example, where the internuclear 
coupling \V^^\ ~ 10~^ K typically, ie., up to 8 orders of magnitude smaller 
than the hyperfine couplings {f|, Fj^}. 

(ii) The Hilbert space of an oscillator is infinite-dimensional and thus quite 
different from that of a spin (or any environmental mode with only a few 
discrete levels). This would not be important if the weak-coupling assumption 
(3) were satisfied; but in the common case where it is not, one has to use spin 
propagators in dealing with the environmental modes. 

Readers working in different fields will recognise many of their favourite mod- 
els in these general forms (quite how general they are has been discussed at 
some length in the literature [12,14-16]). The models are completely specified 
once (a) the Central system Hamiltonian has been given, and (b) the various 
couplings (or more usually, their distributions over q and k) are known. At 
this point two simplifying features usually emerge, viz.; 

(a) The distribution of couplings to the environmental modes shows a fairly 
typical pattern, (see Fig. 1). Coupling to oscillator systems is usually much 
weaker at low energies and temperatures, as the available phase space disap- 
pears, except for the case of "Ohmic baths" such as conducting electrons in 
a metal (and even this coupling rapidly disappears at low temperatures if the 
metal goes superconducting, so that the electronic spectrum becomes gapped). 
Phonons and photons in particular have very weak low frequency couplings to 
most central systems of interest. On the other hand, the coupling to spin baths 
is often much stronger at low energies - this is particularly true for coupling 
to nuclear spin systems, which becomes much stronger as one goes to energies 

Eo = Efc(-^i')^]^^^- For this reason, spin baths are often much more serious 
sources of low-T decoherence than oscillator baths. 

(b) At low temperature the Hamiltonian of the central system is often al- 
most trivial. In microscopic systems like 2-lcvel atoms this may not come as 
a surprise (although an atom is of course already a very complex many-body 
system). In mesoscopic or macroscopic quantum systems the central Hamil- 
tonian is typically describing a collective coordinate which may have a quite 
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large "mass" or "inertia", but whose dynamics simplifies at low energy be- 
cause it is restricted by some confining potential. Thus one ends up with a 
few "canonical models" [17] which have been applied to many different phys- 
ical situations- these include a 1-dimensional harmoic oscillator [18], a simple 
tunneling potential [11,16], and a variety of problems in which the central 
Hamiltonian Ho describes one or more 2-level systems. In fact the problem of 
a single "central spin" (or "qubit") is so important that we pause to sketch a 
few details. 

(iii) Spin-Boson &; Central Spin Models: If we couple a central two-level 

system to either an oscillator or a spin environment, we get two very useful 
models. In the case of an oscillator bath we get the spin-boson Hamiltonian 
[15,16] 

HsB = At, + Cf, + ^[c|f, + (c,^f + + h.c)] + (7) 

with Hos given in Eq. (1). The spin r is in a longitudinal bias field with 
tunneling matrix element A between states | ]) and | |); the couplings c|, ~ 
0{N^^/'^). This model has been widely used [15,16,19] to discuss systems 
like SQUIDs, or nanomagnets, or other 2-state mesoscopic systems coupled 
to delocalised excitations like electrons, phonons, magnons, photons, etc (it 
was originally introduced to discuss problems like nucleons [20] or the Kondo 
problem [21], involving microscopic 2-level systems). 

If we couple to a spin bath we get the central spin Hamiltonian 

Hcs = (Af+exp^S.«-'^'= + H.c) + E^K^fc + ^.p (§) 

k 

with Hgp given in Eq. (5). Here, as noted above, uj\ and ock are often indepen- 
dent of Ng. This model has been used to describe systems like nanomagnets 
[19] or SQUIDs [14,22,23] coupled to nuclear and paramagnetic spins. 

The non-diagonal couplings and in these models are typically smaller 
than the diagonal couplings - for a complete discussion of this point, see Ref. 
[14] . However, one should be careful not to drop them when discussing deco- 
herence (see below). 

These 2 models can be combined (ie., a central spin couples simultaneously 
to spin and oscillator baths- this often happens in reality [14,19]), and they 
can also be generalised to discuss 2 or more spins coupled both amongst each 
other and to an environment - we will use such models in sections 5 and 6. 

(iv) Averaging over the environment: To complete the discussion of these 
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models for open quantum systems one must specify how to average over or 
"integrate out" the environmental modes. This has been extensively reviewed 
[16,14,24] so we simply recall the results here. For both spin and oscillator 
baths the best way is to use path integral methods (lending themselves par- 
ticularly to strong-coupling and tunneling problems). For the oscillator bath 
a common method is to write the proagator for the density matrix as a func- 
tional integral over the "influence functional" J-[q,q'], in the form [16,24]: 

K{Qf,Q'f;Qi,Q'.,t) = 

q{t)=Qf q'(t)=Q'j 

J Vq{T) J Pg'(T)e'(^°l«l-^°l«'])j^[g(T),g'(T)] (9) 

q{0)=Qi q'{0)=Q'i 

where So[q] is the action of the free system (no environment) and J-'[q, q'] is a 
functional over all possible paths q'(r), ^'(t) of the system. J-'[q,q'] is usually 
a rather complex object, and the functional integration over it is not easy 
either- what makes it possible at all is the weak coupling to the oscillators. 
The only cases for which the dynamics of p{t) have been estabhshed so far 
are when the central system is either an oscillator itself, a central spin (the 
spin-boson model) a pair of spins (the PISCES model discussed below), or a 
particle moving either in zero potential or hopping on a lattice, in the presence 
of oscillators. 

For the spin bath a nice simplifying feature appears, originating from the finite 
Hilbert space associated with each spin bath mode. One begins by quantizing 
the bath spins along a z axis, and classifying all the 2^ bath states according 
to their total polarisation M along this axis. Instead of functional integrals 
over all possible paths of the system, one has just ordinary integrals over the 
density matrix ff-°\t) of the free system [14]: 

p(g, g'; t) = i^(0)i°(|/)i^(e, M) (0, y, e, M; Q, Q'; t) (10) 

where the parameters in the free density matrix p'^°\t) now become simple 
functions of (i) a phase variable associated with rotations of bath spins 
induced by transitions of the central system (ii) a similar variable y associated 
with spin bath precession in between central system transitions, and (iii) the 
energy bias e acting on the bath spins, and the polarisation M of the bath. The 
averages A^{(1)), A^{y)i and A^{t,M) are just simple weighted integrations 
over these variables (with weighting functions depending on the couplings to 
the bath and the bath temperature). The replacement of functional integrals 
by ordinary integrals makes the calculation of dynamics for the spin bath 
much simpler than for oscillator baths- in fact one can usually find analytic 
answers. One sees this clearly in the example of the central spin model (8) 
(see ref [14] and below). 
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3 Decoherence for single "qubits": Experimental Implications 

Perhaps the best way to understand a formahsm is to look at some important 
application of it to a real system. Wc tlierefore begin with a very important 
one, to experiments on mesoscopic SQUlDs and magnets- both systems are 
strong candidates for qubits in some future quantum computer. In fact many 
such experiments have been done in the last few years, which we can somewhat 
arbitrarily divide up as follows: 

(i) Experiments claiming observation of coherence, or at least interference, 
between mesoscopic quantum states, either in SQUIDs [25,26] or quantum 
magnetic systems [27-29]. 

(ii) Experiments on incoherent tunneling at the mesoscopic or macroscopic 
scales, also in superconductors [30-32] and magnets [33-37]), done in response 
to theory [11,15,38-42] which indicated that even in the presence of a vast 
number of environmental modes or "microstates" , such tunneling was possible. 

In fact the incoherent tunneling experiments were very successful- and amongst 
other things they indicate that the theoretical treatment of the environment 
works pretty well. Thus, the SQUID tunneling experiments of Clarke et al. [31] 
gave amazingly good agreement with the Caldeira-Leggett [11] predictions for 
electronic "oscillator bath" dissipation in SQUIDs. Recent experiments on 
crystalline ensembles of molecular magnets [37,43] (including direct verifica- 
tion of the role of nuclear spins, by isotopic variation [43]) also agreed quan- 
titatively with the theoretical predictions [19,42,44] for nuclear spin-mediated 
incoherent tunneling in these systems. Thus at least for incoherent tunneling, 
the theory seems to work. 

However, as emphasized long ago in the measurement theory literature, inter- 
ference between 2 "classical" (ie., mesoscopically or macroscopically different) 
states is a much more sensitive test of our understanding of decoherence mech- 
anisms. In the same way it severely tests our understanding of superconductors 
and magnets, and indeed of many-body physics as a whole. The coherence ex- 
periments in SQUIDs and magnets are usually analysed using the spin-boson 
and/or central spin models just given, because in both cases the collective 
coordinate dynamics is governed by a "2-weU" potential at low energies when 
ksT < uJo/27r, where Uo is roughly the same as the gap from the 2 lowest 
states in this potential, to the higher levels. In mesoscopic magnets this "spin 
gap" may easily be 10 K; in SQUIDs it is the Josephson plasma frequency, 
which can be made as large as 1 — 2 K (incidentally, Qq does not decrease 
exponentially with the size of the system, as has sometimes been argued [2]; 
in magnets it may hardly depend on system size at all). 

In this section we discuss the mechanics of decoherence for problems like these. 
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returning at the end of the section to make some brief comments about the 
experiments themselves. Note what is meant by coherence in the case of these 
2-level systems - if, say, we start with the 2-level system in the pure state | t), 
then a purely coherent time evolution thereafter is described by the "free" 
density matrix Pij\t), where (see Fig. 2): 

^ |(l + cos(2Ai)) |sin(2Ai) ^ 



y -|sin(2At) i(l-cos(2At)) j 



(11) 



The job of the theory is to predict the evolution of Pij{t) in the presence of 
the environment - coherence between the states | t) and | |) is then described 
by the ojf- diagonal matrix elements pi2 and p2i- So far the experiments have 
not directly tested any decoherence calculations, and indeed the magnetic 
coherence experiments are controversial- see below. However they soon will, 
and so in this section we discuss what theory says for the decoherence dynamics 
of these off-diagonal matrix elements. 

(i) Spin-Boson Model: We include only the diagonal couphng c| for simplic- 
ity, and assume to be specific that the coupling is Ohmic and weak- formally 
this means that the Caldeira-Leggett [11] spectral function J{uj) = Trauo with 
a -C 1. This model then applies directly to the electronic contribution to deco- 
herence in SQUID tunneling between flux states ±0^ with charging energy 
and a given Q-factor [38]; where a — (160^a;o/-E'c)Q~^- The dynamics for this 
model is fairly well-known [15,16]; in the weak-coupling limit one has standard 
exponental damping of the oscillations, with diagonal and off-diagonal decay 
rates given respectively by [15,45,46]: 

rn(A„o = (A,/£;)^(£;) (12) 
ri2(A„ = ^[2e^(o) + ^He)] (13) 



where the rate function 7(-E) is 

-t{E) = T:aEcoih{E/kBT) (14) 

and E — {^"^ -\- A^)-*^/^ with a renormalised tunneling matrix element A^. A 
somewhat peculiar point is that the size of the off-diagonal matrix elements 
is also reduced, by a factor A^/A ~ (A/Ec)°'^^~°', which explicitly contains 
the high-energy cut-off in the oscillator bath spectrum (Fig. 3); similar effects 
can also be seen for a free particle [47]. 
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We can think of these results as calculations of the rates T{~^ = Tn and 
Tq^'^ = Ti2, with T2 = 2Ti for the symmetric (zero bias) case C = 0. Thus for a 
spin weakly coupled to an oscillator bath, the decoherence rate Ta"^ is simply 
related to the longitudinal rate Tf^. This is partly a reflection of the simplicity 
of the decoherence mechanism here- the longitudinal coupling c| couples to 
a "fluctuating bias fleld" (the oscillators). 

Anyone familiar with the usual arguments for the destruction of off-diagonal 
matrix elements as a solution to the "measurement problem" ought to be 
suspicious of this simple relationship between relaxation and decoherence in 
the spin-boson model. In fact it is too simple, indeed rather exceptional- as a 
general rule one expects no particular connection between them. Results from 
the spin-boson model can thus be misleading, so we now turn to the central 
spin results. 

(ii) Central Spin model: From Fig. 1 we learn that most candidates for a 
"qubit" will begin to feel the deleterious effects of the spin bath at sufficiently 
low T. By itself the coupling a;| causes no decoherence at all- except for a 
rather weak contribution coming from the internal dynamics of the spin bath 
(driven by V^k'). The 2 main sources of decoherence are (i) the precession 
of the spins in the spin bath between transitions of the central spin, driven 
by the term uj^, and (ii) transitions in the spin bath caused by the hipping 
of the central spin (coming from the otk term- the probability for cfk to flip 
during a single central spin transition is IcKfcp when \otk\ <S 1. Both of these 
decoherence mechanisms [42,14] involve the accumulation of a random phase 
by the environment, which when integrated over to obtain the reduced density 
matrix of the central spin, lead to decoherence. The extent of this decoherence 
is quantifled using the parameters 

A = ^EI«^I' (15) 
^ k 

«^ = |ekV4)' (16) 

^ k 

(we assume \ah\, |a;j^/a;|| <S 1). Roughly speaking, these dimensionless num- 
bers tell us the mean phase accumulated by each of these mechanisms whilst 
it operates. Thus we expect strong decoherence if A > 1 and/or /t > 1. To sub- 
stantiate this we give the results for the off-diagonal matrix element pi2{t) of 
the central spin density matrix in each case (for the diagonal matrix elements 
see refs. [42,14]), since these matrix elements quantify decoherence effects: 

(a) Topological Decoherence limit: In this limit the main source of decoherence 
is the transitions in the spin bath caused by a transition in the central sys- 
tem ("co-flip processes"); the time-dependent perturbation on the bath spins. 
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coming from the central spin transition, changes the spin bath state, and each 
such transition gives a random phase (mean value A ) to the bath. One then 
has 

oo 27r 

Pi2(t)=i^(0)pg(0;t)= E (17) 

v=—oo Q 

which is easily evaluated using Pi2{4>]t) = {i/2) sin(A(0)t), where the phase- 
dependent tunneling matrix element is found [42] to be A(0) = A cos(0). Only 
terms with odd winding number u contribute, and one finds that 

P(i2)W = Ee"''^''^^'^'>^2.+i(At) (18) 

In the case A — > (ie., the bath decouples from the central spin) this just goes 
back to the free spin result in (11). If A ~ 0(1) or greater, we have 1^12^*)! ~ 
0{e~^), ie., exponential suppression of the off-diagonal matrix element (cf. 
Fig. 4)! 

(b) Orthogonality Blocking limit: In this limit "co-flip" processes are unimpor- 
tant, but the bath spin precession between each central spin transition now 
adds a random phase to the bath part of the total wave-function- when inte- 
grated over this also leads to decoherencc, in much the same way as with topo- 
logical decoherence. This mechanism becomes particularly important when 
the couplings {^1} are large compared to A, in which case it is important to 
distinguish the different polarisation groups in the spin bath. We then have 

N 

pSW = i°(y)pg(y;^)- E ^M(T)(pg(i, AM(y))) (19) 

M=-N 

where the weighting function Wm{T) tells us the fraction of spin baths, in a 
thermally weighted ensemble, to be found in the M-th polarisation group, the 
phase- weighted matrix element is now Am(?/) = AJm(\/2k^) and the average 
is 

oo 

(pg(i, AM(y))) = Sm,o I dye-yp^^{t, AM(y)) (20) 



This result is rather remarkable- since the fraction of spin baths in an ensemble 
which find themselves in polarisation group M = is ~ (27riV)-i/2, we have 
only a tiny contribution to the off-diagonal matrix elements! The physical 
reason for this result is actually very simple- for coherence to proceed, the 
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state of the nuclear bath must not change between initial and final states. But 
it is easy to see that in the present case this can only happen if M = 0, because 
otherwise the system cannot make resonant transitions (recall that {a;|} > A). 
Consider, eg., the case M = 1. If no bath spins flip during the transition, then 
afterwards the central spin sees M = — 1, ie., the total coupling energy, acting 
here like a longitudinal bias, has changed by roughly 2a;| (for more detail on 
this physics see ref. [19], Appendix B). 

(iii) Experimental Implications; Amongst other things these results demon- 
strate the well-known fact that there is no necessary connection between the 
decay of diagonal matrix elements and those of off-diagonal ones. This may 
seem obvious to those familiar with either measurement theory or NMR. In 
NMR one often sees a huge difference in behaviour between Ti and T2; but 
coherent dynamics exists only over timescales < T2 (off-diagonal matrix el- 
ements) , and not for T2 < t < Ti (even though oscillations in the diagonal 
density matrix elements may still be strong). 

However as far as we are aware this point has not been noted up until now in 
the context of macroscopic coherence experiments. In fact all of these experi- 
ments have either looked at densities of states using thermodynamic measure- 
ments [29], or else at the absorption spectrum of the relevant system (SQUID 
or magnet) as a function of frequency and applied field [22,25,27,28]. We would 
like to emphasize that none of these experiments is then directly demonstrat- 
ing coherence- one can easily envisage a situation in which resonant structure 
in either the thermodynamic density of states or the resonant absorption is 
coming entirely from the diagonal matrix elements. Certainly the SQUID ex- 
periments of van der Wal et al. [22], which look at "ground state tunneling" 
(unlike the excited state tunneling experiments of Friedman et al. [25]), give 
strong evidence that such coherence may soon be demonstrated directly- but 
they do not in themselves constitute such a demonstration. As for the mag- 
netic experiments so far, two of these experiments (both highly controversial 
[48,49]) have been done on large ensembles of disordered nanoparticles, and 
although the specific heat experiment [29] has been done on an ordered set of 
nanomolecules, it does not probe the magnetic dynamics directly. 

However what is certainly true is that there is nothing in the decoherence 
results described so far that rules out macroscopic coherence. Decoherence 
coming from the oscillator bath can easily be made very small for either a 
SQUID (where one can make a -C 1) or an insulating magnet (where phonons 
give very little low-energy decoherence) . The results for spin bath decoherence 
seem more worrying but only if A and/or k are large. In a SQUID the best 
strategy to avoid spin bath decoherence is to make A large (as experimentalists 
are doing), control substrate decoherence, and possibly also to increase the size 
of the SQUID ring. In a magnetic system the best way to reduce decoherence 
effects turns out to be using strong transverse fields, which not only increases 
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A but also radically reduces k; again, the theory indicates that coherence 
should be visible. We do not as yet have any reason to disbelieve this theory- 
in fact it has done rather well in predicting the tunneling characteristics of 
these large systems. Thus the short-term prognosis looks good- it would seem 
reasonable to bet that direct evidence for single mesoscopic sohd-state qubit 
operation will be found in the next year or two. What then? 



4 Coupled Qubits, Quantum Measurements, and all that 

The next step up from single qubits is to couple pairs of them and look at their 
"mutual coherence" or "entanglement" , and how it is destroyed by decoherence 
processes. Ine the same way that one can write down general models for single 
qubits coupled to an environment (the spin-boson and central spin models) 
one can do the same for a pair of them- this is the "PISCES model" , to which 
we come below. However it is interesting to realise that the same model also 
describes a measurement system, itself operating in the quantum regime, and 
set up to probe a single qubit. This equivalence is not reflected in the literature, 
where 2 quite different approaches coexist (somewhat uneasily!). 

(i) Quantum Measurements: The quantum measurement puzzle is perhaps 
the most controversial of all those associated with Quantum Mechanics. The 
first problem is of course to say what it is- this depends partly on one's view 
of "physical reality" , as well as on how one formulates the theory. One simple 
approach begins by noting that that in an "ideal" measurement, the coupled 
"system-apparatus" is supposed to evolve according to 



leaving a "post-measurement" state $j of the apparatus in 1 <-> 1 correlation 
with an initial state (pj of the measured system (obviously we can also have 
\&/ = CjXj^j, with Xj = J2jiOiji(f>i). The puzzle is then simply that for 
the apparatus to function as intended, the states must be macroscopi- 
cally different from each other- meaning that (21) describes a superposition 
of macroscopically different states. 

Answers to this puzzle vary widely. Many formulations of quantum mechanics 
(including the Copenhagen interpretation) make measurements a primitive 
and unanalysable part of the theory, and insist that a measuring apparatus 
must by definition be classical. The superposition ^1'/ is then replaced by a 
mixture, the two being distinguished formally by their density matrices: 



(21) 



j 
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pf^i = \ci\''5ij {mixture) (23) 

Various arguments have been given that, FAPP ("For All Practical Pur- 
poses" [50]), or even in principle, this manoeuvre is justified. Most of these 
start from the idea that the apparatus is extremely "complex", and so has 
many microstates, in a very large Hilbert space, associated with each of the 
"macrostatcs" describing a given outcome of a measurement [2,3]. It is also 
often argued that large amplification and very strong irreversibility (and/or 
chaotic dynamics) are associated with the time evolution of the apparatus dur- 
ing measurement. Hence, it is argued, no conceivable measurement performed 
on "system + apparatus" could distinguish p^^^ from p™^* (it would require 
an operator simultaneously off-diagonal in all the relevant states!). Many of 
the early "decoherence" arguments concerning quantum measurements, going 
back some 50 years, make a similar point [1]. From now on we will refer to this 
set of arguments as the "orthodox argument" for the denial of macroscopic 
quantum superpositions. 

At first glance it appears to be a very powerful argument - after all, the 
macroscopic world does appear to be classical, and the number of microstates 
to which a macroscopic "pointer state" couples is very large indeed. However 
the orthodox argument also suffers from a remarkable vagueness- no attempt 
is made to substantiate it by any detailed quantitative theory, susceptible to 
experimental test on a real measuring system. One is reminded of the 1960's 
conceptual art movement (in which the point was not to create a work of art 
but rather to talk about what it would be like if completed) . The problem with 
the argument is obvious- insofar as it denies the existence and/or observabil- 
ity of mesoscopic or macroscopic superpositions, it is in imminent danger of 
contradicting experiment (see previous section) and certainly confiicts with a 
lot of well-established theory! But what precisely is the flaw in the orthodox 
argument? 

(ii) Coupled qubits: In the Quantum Computation literature [4], on the 
other hand, the above orthodox argument is usually ignored completely (ex- 
cept for the occasional reference to decoherence as "an important problem to 
be resolved" , occasionally accompanied by some rough calculation including a 
phenomenological "noise term"). It is very hard to understand this attitude, 
which entirely ignores the microstates so central to both the orthodox argu- 
ment and to decoherence- after all, the whole point of quantum computation is 
to set up extremely delicate superpositions involving many qubits and a huge 
number of states. So how sensitive is the mutual coherence between qubits to 
the environment? 

Theoretical analysis of this question is simplified by the fact that quantum 
computation algorithms can be set up by coupling pairs of qubits in vari- 
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ous sequences [4]. This allows us to discuss decoherence in the measurement 
process using the same model as we use for decoherence in the simplest quan- 
tum computations, a model involving 2 quantum systems plus the quantum 
environment. Let us now recall this model. 

(iii) Effective Hamiltonian: To address the above questions we consider a 
SQUID magnetometer or a nanomagnctic field sensor, in the quantum regime 
{kT < f2o/27r), set up so its 2 low-energy "apparatus pointer macrostates" 
I t), I J,) couple diagonally to another 2-state system (see Fig. 5). Here we will 
assume the the measured system is also a mesoscopic or macroscopic SQUID 
or nanomagnet; the apparatus pointer states will be described by a spin-1/2 
variable ri, and the system states by a spin r2- Obviously such a set-up can 
also be used as a coupled qubit system if the parameters are adjusted suitably, 
and decoherence is sufficiently small in both of them. 

We may represent the low-energy dynamics of this coupled system by a Hamil- 
tonian which generalises to a pair of spins the spin-boson and central spin 
Hamiltonians given above. For simplicity we will consider here only delo- 
calised environmental excitations (ie., "microstates" ) ; it is easy to gener- 
alise the following discussion, when necessary, to incorporate the coupling 
to a spin bath. Then we get the "PISCES" Hamiltonian [51] HpiscES — 
Ho{ti, T2) + Hint + Hos, with Has given by (1), and 



To be specific in the following discussion we have made some restrictive as- 
sumptions, as follows: 

(i) we assume only longitudinal couplings to the oscillators. A more general 
model would have a coupling fCa/s ^"2 ^2 between the 2 spins Ti and T2, but if 
the systems are mesoscopic, the longitudinal couplings usually dominate the 
low-energy limit [51]. We assume only an oscillator bath- and to be even more 
specific we will assume a T- independent Ohmic coupling to the 2 spins. 

(ii) Wc ignore time- dependence in the parameters- this is because we are 
interested here in environmental decoherence (as opposed to that caused by 
this time-dependence, which is of course crucial to the actual operation of a 
pair of qubits). 

A microscopic derivation of this Hamiltonian for the present situation is ba- 
sically just the derivation of the effective couplings, and this we know how to 
do. For a single isolated SQUID one has a coupling to electronic excitations 
[38], as well as less important couplings to phonons and photons- and one 
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has calculable couplings to a spin bath of nuclear and paramagnetic spins, as 
well as charge defects [23]. In a similar way one can treat the coupling of a 
single quantum nanomagnet to electrons, phonons, and nuclear spins in this 
low-energy hmit [19,52]. When two such systems are coupled one has to do a 
little more work - the couplings to the various baths cause a renormalisation 
of the effective coupling between the two spins - but this can also be done for 
these systems [51]. 

Let us now consider the basis of the "orthodox argument" in the light of all 
this. The microscopic derivations actually do all the work for us - they give 
an explicit representation of the "microstates" associated with each apparatus 
pointer macrostate (these being electronic, photonic, nuclear spin, etc., states, 
coupled explicitly to the pointer states). Needless to say, each macrostate 
is associated with a vast number of microstates - however their number is 
not directly relevant to the problem. What counts is how they couple to the 
macrostate- in particular, how each individual microstate is affected by the 
transition | t) I i)- Here several crucial physical points arise: 

(i) The important oscillator bath states are either gapped, by an energy Eq 
typically much greater than the frequency scale cUo involved in the apparatus 
dynamics (eg., the Josephson plasma frequency for a SQUID); this is the case 
for the electronic excitations in a SQUID or the magnons in a magnet. Any 
remaining gapless excitations (phonons, photons, etc.) have very little spectral 
weight and extremely weak coupling to the apparatus macro-coordinate at 
these frequencies (recall Fig. 1). If the electronic excitations are no^ gapped, 
then their effect is drastic- one gets the famous orthogonality catastrophe 
[53,14] which immediately renders either the SQUID dynamics [15] or the 
nanomagnet dynamics [19] both incoherent and strongly damped. 

(ii) Gapped excitations are only very weakly disturbed by the slow apparatus 
dynamics- this is a typical example of quasi-adiabatic perturbation of the 
environmental microstates, which can be formulated to lowest approximation 
a la Born-Oppenheimer (and more rigourously using the effective Hamiltonian 
approach). We can think of the perturbation of the microstates as a slow 
readjustment to a changing underlying vacuum state (this analogy becomes 
precise when the apparatus state corresponds to that of a quantum soliton, as 
indeed it does both for SQUIDs (a tunneling fluxon) and when the magnetic 
transition corresponds to the tunneling of a magnetic domain wall [41]. 

(iii) The real danger comes from low-lying spin bath states, where this quasi- 
adiabatic argument no longer applies. In the case of nuclear spins we are often 

saved by an "anti-adiabatic" argument- the characteristic energy scale of the 
spin bath states (resulting from the interaction itself) is now much lower than 
ujo, and so again, they are hardly affected by the apparatus transitions (the 
quantitative measure of the effect is provided by the parameters A and/or k, 
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whichever is apphcable - see previous section) . In the case where some of the 
spin bath couphngs lo^ are similar to lo,,. wc have the case of "loose spins" 
[42], where these spins are strongly affected by an apparatus transition- the 
apparatus dynamics again becomes completely incoherent and highly damped. 

From these remarks we see that the main flaw in the orthodox argument is sim- 
ply that it never attempts a proper quantitative discussion of the microstates 
it invokes. However, to get further insight into the physics, and to cover the 
problem of coupled qubits, we must go to the dynamics of the PISCES model. 

(iv) Dynamics: The dynamics of this model in the absence of external biases 
was solved by us a number of years ago [51]. It is controlled by the interaction 
/C, the tunneling matrix elements Ai and A2, the temperature T and three 
friction coefficients cci, 0:2 (describing the Ohmic coupling of each spin to the 
bath), and 0:12 (which describes the effective friction acting on the mutual 
correlations of the spins). Away from the perturbative limit /C <^ A, the 
spins can be either totally locked (if /C ^ T, which results in a single spin- 
boson system with tunnehng matrix element A ~ A1A2//C) or, if /C <^ T, 
still in a correlated phase but strongly disordered by temperature. Coherent 
behaviour of the 2-spin complex is however possible in the "Mutual Coherence 
Phase" (T ^ Aj/aj, /C, where j = 1,2). This phase then allows situations in 
which the combined system-apparatus complex is in a a coherently entangled 
macroscopic state; or in the case of coupled qubits, where a computation can 
be done in time scales much shorter than the decoherence time for this state. 
At this point we note an important feature of this model- decoherence effects 
are massively amplified as soon as /C > A^. 

Let us first use the results to set up a measurement scheme treating both a 
measuring apparatus ta JJ-} and a system : {T; i} in ^ fully quantum 
way and in presence of an environmental bath. To be specific, we consider the 
PISCES model, with bare Hamiltonian 

Ho = AatI + A,r! - Url - K, tItI (26) 

such that /C ^ ^ A^, T represents a strong ferromagnetic interaction and 
A^ » A^; > 1. Without couphng, the apparatus is in state | 1]-) and the 
restriction A^ ^ A^ insures that the apparatus reacts quickly to any changes 
in the system. 

With the apparatus initially in the state | ft") and the coupling turned on at 
some time then a state | fl^t) remains due to bias ^a- On the other hand, 
if at time t the coupled state is | fl^j), fast relaxation to the state | J||) then 
occurs (at a rate ~ (A^/a;o)(-ft'/c^o)^"'^^^) and the combined apparatus- 
system is then essentially frozen. Relaxation to the state | W) takes place at 
the extremely slow rate ~ (A^As//Ca;o)(T/c<;o)^°''*~''^ <^ r']} . In effect, an 
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a. 



ideal" measurement of the form [54] 



^>IT> 



^>IT> 



(27) 



has been performed. The difference in time scales is what makes the measure- 
ment possible. 

Now, let us look at another limit of this model, more appropriate to the 
discussion of decoherence in a pair of coupled qubits. All external biases are 
removed and we assume that (a) each system is identical, and (b) that the 
dissipative couplings a = ai = a2 are very small. If we wish to operate the 
coupled system as a pair of qubits, what then will be the decoherence rates? let 
us consider the specific situation where the 2 qubits are in the pure state | fj) 
at time t = and calculate the diagonal elements Pv^v-iif) with ui = {t, j} and 
^2 — {T; i}- Without going into details (which are complex- the time evolution 
of the density matrix involves three charactristic frequencies and 3 different 
relaxation times [51]) one can make the following general remarks. First, the 
general relaxation of the "ferromagnetic" state, as measured from the time 
behaviour of the combination Pii {t) + p^(t), is controlled by combinations of 
the decay rates rii(2Ar,/C) and ri2(2Ar,/C), using the same notation as we 
used above for the diagonal and off-diagonal matrix elements of the single spin- 
boson problem in Eqs. (12) and (13), and replacing 'y{E) by 7(a//C^ + A^). On 
the other hand, the "dephasing" to the states | fl) and | ||) occurs essentially 
at rates rii(Ar,/C) and ri2(Ar,/C). These are also the decay rates for the 
antiferromagnetic global combination P||(t) + Pi^{t). Thus the principal result 
is that the mutual coherence or "entanglement" correlations are very fragile- 
they are rendered incoherent by a combination of environmental couplings and 
the inter-qubit coupling /C itself, in a way very much like the destruction of 
coherent oscillations for a single qubit once a bias is applied (note how /C plays 
much the same role in the 2-qubit entanglement relaxation as the bias ^ does 
in the single qubit relaxation). 

However we also see that if a is small enough (ie., that F/A <^ 1) then deco- 
herence will not be a problem for the operation of a pair of qubits. One can 
go into more details here, looking at the regions in the parameter space of 
variables /C, T, a, and A where a pair of coupled qubits can operate (the "mu- 
tual coherence regime" defined in Dube and Stamp [51]), but this discussion 
becomes rather technical so we refrain from this here. We also note that any 
realistic analysis of this dynamics must also include the spin bath. 

In any case our main point is clear - by fiddling with the parameters of this 
coupled system we can make it go from a coherent regime to an incoherent 
regime. In the latter the apparatus really acts as a measuring device, and 
in the former as a pair of qubits. We note also that a crucial feature of the 
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apparatus when in the conventional measurement mode is that it is sufficiently 
damped to settle into a single "macrostate" . This is of course is hardly new- 
it is follows from the very definition of a "measurement" that this must be 
so! But this does not affect our point- which is that the same system that 
can behave as a measuring apparatus can also, with only small changes in 
coupling constants or fields, show macroscopic superpositions. One may try 
to argue (like Peres [3,55]) that the word "macroscopic" by definition ought 
to imply that macroscopic superpositions are disallowed (so that coherent 
superpositions of SQUID flux states, or magnetic magnetisation states, must 
by definition be considered microscopic no matter how large the SQUIDs, 
magnets, or flux or magnetisation differences may be). But then the whole 
orthodox argument becomes a mere tautology, and explains nothing- as well 
as forcing a very peculiar notion of the term "macroscopic" upon us (one in 
which a system changes from macroscopic to microscopic by, eg. , tuning of an 
external field!). 

We finish by stressing that these considerations are not academic. Experi- 
mentalists searching for large-scale quantum effects in magnets and supercon- 
ductors must push their SQUID magnetometers or nanomagnetic detectors 
into the quantum regime to attain increased sensitivity. They will eventually 
have to take account of quantum correlations in the functioning of the mea- 
suring apparatus. This promises to lead to all sorts of interesting discussions! 
And of course a world-wide effort is on to build solid-state qubits. Thus the 
decoherence story is only just beginning. 
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Figure Captions 



Fig. 1: Typical form for the dimensionless coupling g{uj,T), with u — T, 
between a mesoscopic system such as a superconductor or a nanomagnet, 
and various environmental modes. In this particular case we assume a 3- 
dimensional system with nuclear spin energy scale ^ 3 x 10^ K, and a 
superconducting transition at T ~ 2 



Fig. 2: Diagonal and Off-diagonal of the density matrix for a "Free spin", 
uncoupled from the environment, and in zero bias field = 0). 



Fig. 3: Diagonal and Off-diagonal elements of the density matrix in the spin- 
boson problem with parameter a = 0.05, T/A = 0.5 and Ec = 50A. The 
value pi2{t oo) = (A/A^) tanh(Ar/2r). The system is again unbiaised, ie.. 



Fig. 4: Diagonal and Off-diagonal elements of the density matrix in the central 
spin problem with zero bias, and with topological decoherence dominating. 
The value of the topological decoherence parameter is relatively weak: A = 1/8. 



Fig. 5: The general problem involved in dealing with a measuring system 
(collective coordinate Q' coupled to some other quantum system (collective 
coordinate Q), with both coupled to a quantum environment; all variables are 
quantised. The same kind of model describes a pair of coupled qubits. 
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Fig. 2. 
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Fig. 4. 
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